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q ", Abstract 

Configurations of pure SU(2) gauge field theory on the lattice are transformed 
fNJ . to Landau gauge. After Fourier transformation, large momentum amplitudes are 

\ suppressed by a variable amount, and the configurations are transformed back to 

x-space. Spectacular peaks in electric and magnetic field strengths are found, which 
share many properties with either almost pointlike instantons in regular gauge or 
with extended anti-instantons in singular gauge. Environments around those peaks 
. are visualized with respect to the action, to the topological charge density, to the 

gauge fields and to electric and magnetic field strengths. The density of the peaks is 
of order 1/m -4 , and it scales according to the string tension under a variation of the 
coupling constant. The evolution of the peaks under the amount of Fourier filtering 
is visually compared to the evolution under cooling, and the gauge dependence of 
the peaks is discussed. Advantages and shortcomings of this method are discussed, 
with emphasis on possible strong distortions of the action and the topological charge 
density, which become gauge dependent. Finally, I compare the character of the 
i ' SU(2)-configurations to those of non-compact abelian theory. 
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1 Introduction 

It is hardly necessary to justify an attempt to look directly onto gauge field 
configurations which have been generated by standard lattice techniques. Un- 
fortunately, the amount of difficulty encountered in the course of such attempts 
is also quite obvious. On the level of physics, the freedom associated with local 
gauge transformations requires a careful discussion of the significance of any 
observations, other than those based on a few gauge invariant quantities. Even 
for a given gauging algorithm, the resulting configuration will depend on the ap- 
plication of a random gauge transformation prior to gauging, due to the Gribov 
ambiguities Furthermore, the short distance quantum noise degrades any 
visual impression. The removal of theses fluctuations will, in any case, distort 
the lattice configuration to an unknown degree. 

On the technical level, the number of local degrees of freedom in non-abelian 
theories is so large that one has to strongly reduce the amount of information 
before something can be plotted on a screen or on paper. On the other hand, 
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the speed of progress in computer-based visualization techniques is so high that 
today's obstacles may be irrelevant next year. 

The present approach to dealing with the physics problems is the following: 
A thermalized gauge field configuration is transformed to the Landau gauge, 
and the gauge fields of the resulting configuration are smoothed by standard 
Fourier filtering procedures. Thus, the short distance fluctuations are elimi- 
nated. Presumably, these are of perturbative origin and do not deserve too 
much attention. To what extent this is correct, may be partly checked a pos- 
teriori. There are good indications that the assumption is not totally correct. 
These indications are connected with the observation that at a few positions on 
the lattice, gauging leads to very strong fluctuations of the gauge fields which 
look like singularities, and that at these points the Fourier amplitudes at all 
momenta are non-perturbatively connected. This important observation will be 
explained in detail in sect. [42], where also the nature of these special points 
will be described thoroughly. 

In order to avoid being overwhelmed by the amount of data on a large lattice, 
a subset of the lattice points is selected and the associated fields are visualized 
on a PC-screen. The selection of what to plot and what to discard is partly 
made dynamically, and the techniques described in the following are motivated 
by the study of instantons on the lattice. After only quite modest filtering, 
i.e. after discarding just the highest momentum Fourier components, one can 
observe pronounced maxima in several quantities such as electric and magnetic 
fields, in the Wilson action density^] S a (x) and in the topological charge density 
P a (x). Their positions are identical with those of the gauge field singularities 
mentioned previously. They are found numerically by a computer sweep through 
the lattice. After this, the selection of lattice sites is done either by considering 
three 2-dimensional planes through one of those maxima, or by selecting a 4- 
dimensional subcube (e.g. of size 9 4 lattice unites) around the peaks. 

The density of the peaks agrees with that of instantons found in the SU(2)- 
vacuum by cooling f2|, and the visualization tools described below suggest that 
we are indeed seeing instantons on the lattice. Especially, the orientation and 
the relative strengths of electric and magnetic fields fit well with those of (anti-) 
instantons. However, using the action density S(x) as an indicator, the size 
of the objects seems to be significantly smaller than the size of the instantons 
obtained by cooling |§. In this connection it has to be realized that the effect 
of Fourier filtering on gauge configurations will be quite subtle when the gauge 
fields are so large or so slowly varying that the non-linear properties of field 



strengths are crucial. This will be discussed in sect. ET2L where the influence 



of Fourier filtering on artificial lattice instantons is also studied. This influence 
turns out to be very strong if the instanton is realized in the singular gauge (see 
ref. || and eq. |20] for a description of this gauge). Now, it is quite difficult to 



1 The upper index a refers to the 3 colour components of SU(2). 
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decide unambiguously -on the basis of the lattice configuration- whether we are 
dealing with an instanton with a radius of one or two lattice units in the regu- 
lar gauge, or with a broad anti-instanton in the singular gauge. This difficulty 
is connected with deviations of the observed gauge fields from the ideal field 
distribution. Thus, the visualization technique provides us primarily with reli- 
able information on the rapid variation of the gauge fields at the special lattice 
points, whereas an identification of the fields with the pure case distributions 
will require more detailed work. 

Returning to the description of technical tools, I list the following methods: 

• A projection of all local maxima of e.g. the topological charge density onto 
one plane, where the maxima are defined within certain small hypercubes. 
A first impression is given by fig. [I], derived from a configuration on a 
48 3 x 64 lattice at f3= 2.85 which shows the presence of 5 clean peaks and 
one murky one with (anti-) instanton properties, plus some unqualified 
noise. This plot is derived via a relatively weak Fourier filtering (F.f.). For 
the effect of stronger filtering, see fig. The details of this pictures are 
described in section £|. 

• Gauge fields, electric and magnetic fields are plotted in three lattice planes, 
passing through one of those peaks. 

• Field lines of electric and magnetic fields passing through the maxima are 
plotted for the three SU(2) colour components simultaneously. This allows 
one to study the quality of the wave function of the instantons as function 
of the distance from the center. An example can be found in fig. ^|, where 
electric fields are plotted by 3D techniques. Rotation of the picture on the 
screen shows that the field lines of the three colours are orthogonal on each 
other, as it must be for an instanton. 

• The evolution the gauge fields on a small sublattice, as a function of the cut- 
off, inherent in the F.f., can be visualized on the PC-screen (not on paper) 
by means of animation techniques. This allows to demonstrate how the 
maxima emerge out of the quantum noise when the cut-off in momentum 
space increases beyond a small fraction of an inverse lattice unit. The 
same animation can be performed for cooling, which demonstrates the 
qualitatively different response of the configuration on the two techniques. 
The outcome is described in sect. [5]. 

Further investigations include the following: 

Fourier transformation requires some kind of linearization. A comparison of 



two different linearization techniques which are described in sect. [2.1| , indicates 
that the phenomena described below do not depend on these technical details. 

Since the gauging procedure is not unique |Lj , an ensemble of Landau gauges 
is studied by first performing several random gauges on a given configuration. 
This gives some insight as to what depends on the special Gribov copy into 
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which the gauging has led, and what not (see sect. ^). A comparison with the 
results of other smoothing techniques, notably the cooling method[0, || 0], is 
advisable and will be attempted to a certain degree in sect. 

The interpretation of individual plots is, of course, open to subjective judge- 
ments. It has to be supplemented by measurements with sufficiently good statis- 
tics, before definite conclusions can be drawn. This criterion is perhaps not fully 
met in this paper, and the work has to be regarded as exploratory. 

Technical information on the visualization and on the access to coloured 
figures is given at the end of section [| 

2 Fourier Filtering of Gauged Configurations 

It is well known that lattice configurations are so rough that gauge invariant 
quantities like the Wilson action vary rapidly on the scale of one lattice unit. 
This has already been visualized by refs. Q and for the action and the 
topological density, and it holds even at the largest (3 accessible to sensible Monte 
Carlo simulation at present. It is natural to expect that in a smooth gauge - 
like the Landau gauge- the gauge fields themselves share the above property. 
Thus, before anything can be seen in a lattice configuration, a modification of 
the configuration is necessary, and the ensuing distortions may be quite severe. 
In this paper I concentrate on Fourier filtering (F.f.) as a method of removing 
short range quantum fluctuation. The filtering will act on gauge fields which 
are obtained by linearization of lattice link variables in the Landau gauge. In 
this gauge, all link variables are transformed towards the continuum limit as 
closely as possible, with no position or direction preferred. An extension of 
the visualization techniques described here to non-covariant gauges is certainly 
possible, and it can be pursued at a later date. 

It is not being claimed that this method is really superior to other smooth- 
ing techniques. However, it is conceptually quite different from the standard 
cooling technique [0, R M, [J, and, what's more, it leads to different results. It 
is suggested, therefore, that it may be advantageous to have access to different 
types of tools when discussing in what sense something exists on the lattice. 

2.1 Linearization and Selection of Variables 

Before discussing possible advantages of the F.f., its implementation in a non- 
linear theory has to be defined. Naturally, in the course of this, severe ambi- 
guities will become apparent. I will work with the representation of the SU(2) 
link variables 

U^(x) = u 0>ll (x) + itu^x) (1) 

with 

< M (z)+^(x) = l. (2) 
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Before being able to use a Fourier transform, I have to define some linearization. 
I first tranform to Landau gauge by numerically searching for a maximum of 
the gauge function F(U), 

F(U) = J2u ,,(x) (3) 

under local gauge transformations. In the weak coupling limit, all "large" com- 
ponents uo^(x) should approach unity, which would easily allow linearization 
by working with the "small" components only. In practice, however, even at the 
largest values of (3 used here, many link variables are far from this weak coup- 
ling limit. Especially, at f3 — 2.7 with the Wilson action, one finds a fraction 
of 5 x 1CT 6 links with uo tll (x) < 0. This poses a difficulty for the linearization 
which is necessary for a F.f. 

I have used 2 different approximations for dealing with this departure from 
the weak coupling regime. The first one is to simply identify the (non-compact) 
gauge fields with the "small" components of U^x), i.e. 

2 

-u^x) -> A^x). (4) 

— * 

After modifications of the A fJi (x) by F.f., link variables have been reconstructed 
by setting u'^x) = ^A' Ax) and 



u'o,^x) = ^Jl-u^x), (5) 

if u 2 (x) < 1. Otherwise, I set u' (x) = and normalize the u' (x) accordingly. 

Obviously, this is a poor approximation to the link variables if u ^ < since 
an adiabatic approach to the continuum limit u jM 1 should proceed first via 
values uq^ ~ and thus via increasing values of u 2 ^(x). Normally, a F.f. will 
proceed via decreasing amplitudes almost everywhere. 

In order to obtain a better mapping for the region uq ~ 0, I use a kind of 
stereographic projection 

A, = 1^2/(1 + ^)- (6) 
9o 

After an eventual smoothing modification of the A^ by F.f., normalized link 
variables can be reconstructed via the obvious inversion: 

* = %X>, (7) 



u' = Jl-w 2 /4w, (8) 



u' = ±\]l-u\ (9) 

where the positive root in the last equation is taken for w 2 < 2, and the negative 
one otherwise. This transformation brings the "south-pole" u° = — 1 not to 
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infinity, but to w 2 = 4. In section 4.1 I will describe the minor differences due 



to the two techniques. It should be stressed again that the prescriptions are 
arbitrary, similar to all other smoothing techniques. 

After having defined the gauge potentials A^, lattice Fourier transforma- 
tion leads to momentum space amplitudes A^(k), and these are modified under 
filtering by the substitution 

A'/k) = exp(-X 2 k 2 )A,(k), (10) 

where the filtering parameter A 2 will be varied in a large region which is deter- 
mined by the momenta available on the lattice. These are, for twisted boundary 
conditions in x, y and z directions (with k = 2 sin |): 

2vr(n + i) , 
k x , n = T - , n = 0,...,L x -l (11) 

for colour index 2 and 3 (or green and blue), and for colour index 1 (or red) one 
has 

2tTTI 

k x , n =-j—, n — 0, L x — 1. (12) 

Corresponding expressions hold for other directions. Thus, the minimal squared 
momentum is not zero but, for the large lattice of size 48 3 x 64, = 2n 2 /L 2 . = 
0.0086/a 2 . In this value of A 2 = 10 a 2 will eliminate almost all momentum 

space amplitudes, and such a cut-off is the largest one which will maintain at 
least some residual structures. 

After transforming the fields back to x-space, electric and magnetic fields in 
direction i are defined by open plaquettes: 

Ei{x) = -U e>i (13) 
9o 

with 

U e4 (x) = U t (x)Ui(x + e t )Uj(x + ei )U}(x), (14) 

and a similar expression for the magnetic fields Bi(x). 

Variables to be plotted are thus the coloured v4 M (a;)-fields, the electric and 
magnetic fields and x-space scalar products for the three colours individually: 
The three components of the action density, 

s a (x) = i(E (^ a W) 2 + E (B-(x)) 2 ) (15) 

Z 8=1,3 8=1,3 

with a= red, green and blue, and the three components of the pseudoscalar 
density, 

P a (x) = £ E?Bf(x), (16) 

j=1.3 
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which are related to the topological charge density q(x) by 

= j| E pa W- (17) 

a=r,g,b 

The topological charge is 

<? = $>(*), (18) 

a; 

which is integer-valued in the presence of clean (anti-)instantons. Visual results 
and semi-quantitative results for these variables will be presented in section [|. 



2.2 The Pro's and Con's of Fourier Filtering 

Ideally, a smoothing technique like F.f. should make significant structures in a 
lattice configuration somewhat broader, while reducing their amplitude more or 
less monotonically as a function of the cut-off parameter A 2 . At the same time, 
it should decrease the amplitudes of insignificant fluctuations to a large extent. 
This expectation may work out for the gauge fields A^x), but not necessarily 
for derived quantities like the field strengths and the action. The result of 
smoothing will also, under certain conditions, depend sensitively on the gauge 
even for quantities which are gauge invariant before smoothing. This holds true 
because the concept of momentum space gauge fields clearly is gauge dependent. 
Possible effects will be discussed in connection with instantons since instanton- 
like configurations turn out to be the most impressive objects for visualization. 

Important effects of F.f. have two sources. First of all, the field strengths 
in non-abelian gauge theory are, at certain distances from the center of an 
instanton, highly sensitive to a cancellation between the abelian term d^A^(x) — 
d u A^(x) and the nonlinear term gQ€ab c A b (x)A^(x). In the regular gauge, this 
cancellation occurs at large distances from the center. The effect of F.f. then is 
that there the action increases for increasing A 2 , whereas at small distances the 
opposite holds. 

Secondly, the gauge fields may be so large that any linearization technique 
will provide us with misleading results for the field strenghts and for the action. 
This will happen e.g. for instantons in a singular gaugef] or for almost pointlike 
instantons. Both the abelian term and the nonlinear term of the field strengths 
may then be distorted strongly by F.f., such that the results will be rather 
surprising. 

Since the magnitude of such effects will depend on the instanton radius p, 
one has to specify p for a quantitative discussion. According to ref. 0, we 

2 On the lattice, such instantons will be defined by converting the gauge fields of a regular instanton 
to lattice variables using the stereographic method eq. ^, by transforming the gauge fields of the singular 
gauge transformation (see eq. [2^ below to lattice variables in the same way. Finally, this (almost singular) 
gauge transformation is applied to the regular gauge distribution. 
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expect to see around 1.5 (anti-) instantons per fm A with an average radius of 
< p >~ 0.2fm, which translates to < p >= 4.5 a at (3— 2.7 and to < p >= 7.1 a 
at (3= 2.85. Thus, it is interesting to study the effects of F.f. on an instanton - 
which has been artificially set on the lattice - with p = 5 a. 

In the regular gauge, the effect of a cut-off A 2 = 0.5 a 2 leads to a reduction 
of the action density at the center of the instanton by a factor 1.4. At large 
distances, the action increases slightly upon filtering which is contrary to expec- 
tation and which signals a possible weakness of the method. However, since the 
action is small in this region, the effect is not spectacular. Close to the center, 
the cancellation mentioned above is insignificant because the gauge fields are 
small, but vary rapidly. Thus, the abelian term dominates, and there the result 
of smoothing is as expected. Finally, we note that the distribution of the gauge 
fields changes quite moderately. 

However, if the (artificial) instanton is realized in the singular gauge, the 
results after smoothing for the action and for the pseudoscalar density may 
be quite spectacular. Before demonstrating this, I quote the singular gauge 
transformation ||. The standard instanton gauge fields in the regular gauge 
are 



with the 't Hooft coefficients [|K| rj afiv and 77 (used below). The singular 
transformation 



Vr 

leads to gauge fields which fall off more rapidly at infinity, but which are singular 
at the origin: 



Of course, on the lattice there is neither a real singularity at short distances, 
nor at infinity. 

Now, under the conditions mentioned above, the action develops a sharp 
peak, and the sign of the pseudoscalar density within the region of the peak is 
reversed after F.f. The rapid fluctuations of the singular gauge fields, around 
the center of the instanton, get weakened by F.f., but they are still clearly visible 



The sign flip can be understood easily. The gauge field distribution of a broad 
instanton in the singular gauge is not too different from the distribution of a 
regular anti-instanton distribution of radius p ~ la, as long as only a restricted 
interval of distances is considered. It is the weakening of the short range spike 
in the gauge fields by F.f. which converts the singular gauge instanton into a 




(19) 




(20) 




(21) 



for A 2 = 0.5 a 2 . 
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regular anti-instanton - with opposite pseudoscalar density - for a limited range 
of distances. 

Another shortcoming of gauging is the Gribov problem |l|]. On the lattice, 
it manifests itself in the fact that any gauging algorithm which searches for the 
maximum of F(U), see eq. ||, will end up in one of many local maxima, at 
least on a large lattice. The gauge fields at these local maxima are distinct, and 
this affects an eventual view on the configuration, as well as average quantities 
like the gluon propagator. In this investigation, it has been checked to what 
extent certain structures are strongly gauge dependent in the sense that they 
appear only in one special local maximum and not in others. Before doing so, 
the structures will have to be investigated. They will be described in section f|, 
and the test is contained in section 

On the positive side of F.f., the separation into low and high momentum am- 
plitudes seems to have principle advantages. Asymptotic freedom in non-abelian 
gauge theories implies that the momentum amplitudes close to the lattice cut- 
off can be treated perturbatively, i.e. they should be independent of possible 
non-perturbative configurations at low momenta. It is not clear to what extent 
this is correct. In section [|]2] it will be shown that the excitation of amplitudes 
with momenta of O(10 GEV) leads to a clustering of gauge fields around the 
positions of non-perturbative objects, presumably of the instantons. Thus, a 
possible failure of the above assertion can be studied by F.f., and this should be 
regarded more as an advantage than as a drawback of the method. 

Other advantages of the Fourier filtering technique have to be discussed in 
comparison especially to the well known cooling technique, and in the investiga- 
tion of special objects. Cooling drives the configuration to local minima of the 
action, and thus it may create objects which "exist" on the lattice for a certain 
number of cooling steps and disappear again in the process of further cooling. In 
F.f., such objects can be attributed to low momentum components. They may 
be masked by high frequency components and become visible only after removal 
of a certain amount of high momentum amplidudes, but on further smoothing 
they just get broader and will not disappear. This will be demonstrated in 
section [|. 

It has to be noticed that cooling induces a more violent distortion of the lat- 
tice configuration than F.f. By applying 0(30) cooling steps to a configuration 
after gauging, the low momentum amplitudes are suppressed significant lyQ, i.e. 
by 20% or 30%. In F.f., nonperturbative objects become visible for a cut-off 
A 2 ~ 0.2 a 2 , i.e. if momentum amplitudes with k < 0.5 a" 2 are modified by less 
than 10%. For such momenta, the suppression of amplitudes by cooling already 
amounts to a factor 3 or more. 

Finally, the difference of F.f. as compared to other methods is interesting in 
itself. It may be that the comparison of various techniques helps to clarify in 

3 Low momentum amplitudes are again denned by Fourier decomposition. 
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which sense non-perturbative objects exist on the lattice. 



3 Visualization Techniques 

The lattice variables described in the previous section are too numerous to 
allow a complete presentation on a PC-screen of present quality. This holds 
for the set of lattice points in x-space, for colour and for the Lorentz index. 
Therefore, I select a lattice point which is suspected to be interesting according 
to some criterium, predominantly because the Wilson action density has a local 
maximum (after F.f.) at this point. Then, certain environments around this 
point will be selected. These may be 

• three 2D-planes (always including the t-direction) through the special point, 

• a small (e.g. 9 4 lattice units) subcube 

• electric or magnetic field lines, emanating from the special point. 

Another technique is to select subcubes around several (e.g. 50) maxima and 
project the variables into lower dimensions. Details will be given below. 

Once the variables have been selected and written to some file, a FORTRAN 
routine reads this file, calculates 3D-angles and colours and writes a VRML 
("Virtual Reality Modeling Language", see ]TT[) file. In detail, each oriented 
variable is represented by a coloured 3D-cone, pointing into some 3D-direction 
which will be specified by the variable. For S a (x) and P a (x), this direction 
is given by the orientation in colour space, whereas for electric and magnetic 
fields it is given by the direction in x-space. In many cases, the SU(2)-colours 
are mapped onto the ordinary PC-colours red, green and blue, but colour can 
also be used to indicate positions in the 4th dimension. For the gauge fields 
Afj,(x), the situation is not so straighforward. Whereas the SU(2)-colours will 
be represented as before, one of the 4 space-time directions has to be suppressed. 
Which ones are selected seems to be rather arbitrary, and various choices have 
been tried. 

The size of the cones is given by 

• the sum over colour indices in the case of S a (x) and P a (x), 

• or by the spatial length for Ef(x) and Bf(x) (SU(2)-colour a fixed), 

• or by the 4D- length in the case of A a ^(x). 

The VRML-file can be interpreted and sent to the screen by a standard 
browser, e.g. by the well known and free Cosmo-Player [F2|. I use the freely 



available browser GLVIEW [|13[] with good results. Both tools work for Win- 



dows NT, and the search for a LINUX-browser is continuinefl GLVIEW allows 



4 There seems to be much activity in the internet in connection with a LINUX browser, but nothing 
definite and easily implementable is known to me. 
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to shift, zoom and rotate the ensemble of cones, and thus it enables a closer 
inspection of interesting patterns. The capacity of the browser is sufficient to 
plot a complete plane of a 48 3 x 64-lattice with 3 SU(2)-colours, e.g. for the 
electric fields. In this way, one easily recognizes the correlation of maxima in 
field strength between the three colours, as they are typical for instantons. 

Another important feature of VRML is animation. This allows to visualize 
the evolution e.g. of P a (x) as a function of increasing A 2 . For this, around 
10 different values of A 2 between A 2 = and A 2 = 0.5 a 2 are used to generate a 
sequence of small planes around a previously detected maximum of P a (x). 

True 3D representations will certainly be feasible with the help of special 
hardware, but this will not be pursued here. 

The coloured figures are enclosed in the format .png, and the figures 4 and 
18, which are black/white, are enclosed as .ps-files. The reader is also referred to 
my homepage www.desy.de/~gutbrod, where, in the link "publications/visuali- 
zation, html", one can study and print the figures in the .jpg format, albeit in 
somewhat lower quality. 

4 Results of Visualization 
4.1 Global Plots 

I will start the presentation of selected examples by an overview over full lat- 
tice configurations, where local maxima of the action density have been used 
to identify "hot spots" . In detail, the lattice of size 48 3 x 64 is subdivided into 
sublattices of size 8 4 , and positions of the local maxima of the Wilson action are 
determined within these sublattices. This is done after Fourier filtering with a 
cut-off in the range A 2 = 0(la 2 ). If a pair of maxima has a distance less than 
3 a, the lower maximum is discarded. The lattice positions of an environment 
around the maxima are stored, as well as the corresponding values of the pseu- 
doscalar density P a (x) (without colour summation). In fig. [I], P a (x) is plotted, 
with 3 space-time dimensions properly respected. The fourth dimension, z, 
is converted to colours interpolating between red, corresponding to z = 1, and 
green, corresponding to z = L z = 48. The y— direction becomes visible properly 
only when the plot is rotated around the (vertical) X cLXIS. This is accomplished 
with the help of the VRML-browser. The direction of the cones on the plot is 
given by the direction of P a (x) in colour space. 

The configuration belongs to an ensemble of 10 thermalized configurations 
at (3 = 2.85, with twisted boundary conditions 0. The cut-off parameter A 2 
has been chosen as A 2 = 0.3 a 2 . In fig. ^| the same quantities are plotted, 
but with A 2 = 1.0 a 2 . The overall character of the plots is the same for all 
configurations. Intensive investigations with different lattices sizes, with values 
of (3 in the range 2.5 < (3 < 2.85 and with random gauge transformations before 
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the Landau gauging procedure have shown that the character of these plots is 
by no means accidental. 

The surprising structure of the pictures, especially that of fig. |], is the clear 
separation between six pronounced spikes and a noisy background, without 
any objects with interpolating character. Four spikes are nicely parallel to the 
diagonal in colour space, and two are antiparallel. Within the spikes, the action 
density equals the pseudoscalar density within 20% (this is not shown). More 
detailed investigations, which will be explained in the next section, show that 
the spikes share many properties of (anti-) instantons. Surprising facts are 
their small radii and an early corruption of the instanton wave functions as a 
function of distance from the center of the peak. Anticipating the outcome of the 
discussions given below, it is more likely that the objects are broad instantons 
in a singular gauge than almost pointlike instantons in a regular gauge. 

Only minor details of the plots depend upon which of the two methods of 
linearization, described in section pTTL was used. Especially within the spikes, 
the pseudoscalar density may differ by 10% between the two methods, and by 
much less outside the spikes. 

The suppression of a transition region between the spikes and the background 
can be made more quantitative, and, in connection with this, the density of 
the spikes and its scaling behaviour can be determined. On each of the ten 
configurations, I have identified 0(1.500) local maxima of the action density. A 
value of A 2 = 0.5 a 2 has been used. Then, the action has been summed within a 
radius of R = 2a around the positions of the maxima which are called x m . The 
resulting contribution to the action will be denoted by 

(x-x m ) 2 <R 2 

£ m (A 2 ,i?) = £ S ( x )- ( 22 ) 

X 

The number of £ m (A 2 , _R)-values found in intervals of unit length, summed over 
all 10 configurations, is shown in fig. | as a histogram with full lines. Because 
of the clear separation of the peak in the histogram at high values from the 
background, the number of spikes per configuration has a relatively small sys- 
tematic error. One finds 6.1 ±0.8 spikes per configuration, which transforms to 
a density in physical units 

D sp {(3 = 2.85) = 1.4 ±0.3/m" 4 (23) 



if the value of a((3 = 2.85) = 0.028/m of ref. |TJ]] is used. This density is 
consistent with the number of instantons given in || at smaller values of (5, 
which has been obtained by the cooling method^. A similar density has been 
found for SU(3) in ref. ||, but it is smaller by a factor five than the result for 



5 This agreement does not extend to the observed radii. As will be discussed below, the spikes seen in 
F.f. have a radius of one or two lattice units, whereas in || an average radius of 8 lattice units has been 
found. 



12 



SU(3) obtained in by somewhat stronger cooling. It has to be mentioned 
that predictions for the instanton density play an important role in the search 
for non-perturbative phenomena in deep inelastic scattering ]15|. . Before the 
present result with its relatively small error can be used in this context, a better 
understanding of the nature of the peaks is mandatory. 

Also included in fig. [| is the summed action S(A 2 ,i?) for 20 configurations 
on a lattice of size 32 3 x 64 at (3 = 2.7. This histogram is marked by a broken 
line. The corresponding plot of one configuration with spikes is shown in fig. [5| 
It is obvious that the character of this configuration is slightly different from 
the one shown in fig. At the lower value of (3, a transition region between 
pronounced peaks and the background is clearly visible, and the same can be 
observed in the histogram. 

For a scaling test of the density of peaks, the vertical scale in the histogram 
for f3 = 2.7 has been scaled by the ratios of the volume (in physical units) in 
comparison to the case (3— 2.85, i.e. by 

V phys ((3 = 2.85) = V latt ((3 = 2.85) a 4 ((3 = 2.7) 

V phys ((3 = 2.7) V latt (/3 = 2.7) a*((3 = 2.85)' 1 ' 

The size of the lattice constant a((3 = 2.7) has been taken from ref. [|l7j]. This 



leads approximately to the same area under the spike region. More quantita- 
tively, at /3 — 2.7 there are 13.4 ± 2 spikes in a configuration, and we find 

D sp (f3 = 2.7) = 1.6 ± 0.4/m" 4 . (25) 

If one derives the ratio of the scales from the density of spikes per configuration 
at (3 = 2.7 and (3 = 2.85, one obtains 

a(2.7)/o(2.85) = 1.64 ± 0.08, (26) 



which is in very good agreement with the result from the string tension p4| , i.e. 
a(2.7)/a(2.85) = 1.60 ±0.06. 

Scaling of the density alone does not necessarily imply that we are dealing 
with physical objects. If the spikes were instantons in a regular gauge with a 
radius of a couple of lattice units (assume, for definitenes, a radius p = 5 a at (3= 
2.7, corresponding to a value p = 0.22 fm), then at (3= 2.85 the radius should 
have increased at least by a factor 1.5, i.e. to 7.5 a. Since we sum the action only 
within a radius of 2 a, the quantity S(A 2 , R) would essentially be proportional to 
the action density of the instanton around the center. This, however, decreases 
as p 4 / (x 2 + p 2 ) 4 , i.e. rather strongly with p for x 2 <C p 2 . For a scaling behaviour 
of the radius, the peak of the histogram should have been shifted to the left by 
a factor 4, when going from (3= 2.7 to (3= 2.85. Within our statistics, the peak 
has moved to the right. This excludes scaling of the radius, and there is good 
evidence that the objects have a pointlike character. It has to be remembered, 
however (see section |2.2| ), that extended artificial instantons in a singular gauge 
behave abnormally under F.f., as the action develops an unphysical spike, and, 
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within this spike, the pseudoscalar density reverses sign. At present, I see no 
convincing method to resolve this ambiguity. A safe conclusion is that in the 
Landau gauge the short range spatial variations of the gauge fields cannot be 
described fully by perturbation theory. Some lattice points, with a density of 
~ 1.4//m 4 , show a violent oscillation even after modest F.f., which leads to 
significant peaks in the field strengths, presumably via the term d^A^ — d v A a ^ 
without compensation by the non-linear term. The detailed properties of these 
variations will be studied in the next subsection. 

4.2 Properties of the Spikes 

The visualization technique can be used to study how well the spikes, observed 
in the global plots, share properties of instantons. These properties are 

1. a strong localization of the action density in space-time, 

2. selfduality, i.e. = ±B?, a — r,g,b, i = x, y, z, 

3. colour democracy, i.e. all colours are excited with equal strength, and 
electric (magnetic) fields of different colours are mutual orthogonal, 

4. and a special dependence of the gauge fields on the distance from the 
center, which leads to a cancellation between the linear and the non-linear 
terms in the electric and magnetic fields. In the regular gauge, this occours 
at large distances. 

In fig. I show the pseudoscalar density on a lattice plane running through 
one of the spikes of fig. 0, with a cut-off A 2 = 0.15 a 2 . It is impressive to note 
the strong localization and the vastly different level of the spike amplitude, as 
compared to the background intensity. It is also obvious that F.f., with the 
relatively low value of the cut-off, does not really smoothen the background, 
which is still quite noisy. Primarily, the F.f. reduces the background intensity 
such that it does not mask the peak. In addition to the strong localization, 
the plot also shows that the direction of the spike is parallel to the diagonal in 
colour space. This is consistent with colour democracy. 

In fig. |7|, the gauge field distribution in the same plane as above is shown for 
the red colour, where the fourth space-time component is disregarded in calcu- 
lating the direction of the conesf]. At the position of the spike, one notices a 
slight increase in the activity of the gauge fields, which seems to be statistically 
barely significant. But, checking for the other colours and for more configura- 
tions, the same phenomenon always shows up, and the correlation between the 
magnitude of gauge fields and the spike in the action has to be accepted as real. 

The observation of a peak in the gauge fields allows to ask how well compo- 
nents of high momentum are determined by perturbation theory. For a cut-off 
A 2 = 0.5 a 2 , the peak in the gauge field intensity is reduced by about a factor 

The size of the cones, however, is proportional to (XL=o z(A r ^ d )' 2 ) 1 /' 2 . 
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2. This has important consequences. It implies, that those Fourier components 
which get removed when proceeding from A 2 = 0.15 a 2 to A 2 = 0.5 a 2 , are cou- 
pled dynamically to the non-perturbative object which is responsible for the 
spike. In momentum space, these components have a momentum k 2 « 2/a 2 ~ 
120GeV^ 2 /c 2 . Obviously, as already mentionend in section |2.2] , Fourier ampli- 
tudes around these momenta cannot be determined completely by perturbation 
theory. 

The quality of the wave function of the hypothetical instantons can be stud- 
ied by plotting e.g. electric field strengths. In fig. |8|, I show the coloured 
electric field strengths, attached to lattice points in the same plane as above, 
but with a smaller environment around the spike selected. The cut-off is again 
A 2 = 0.50 a 2 . One nicely observes how the maxima of absolute values pop up in 
the same region for all colours, and how the field strengths of different colours 
are mutually orthogonal. It is evident, however, that there are several other 
locations where the field strengths are non-negligible. The amplitudes of these 
excitations amount to 1/2 or 1/3 of the maximal values at the spikes. The 
nature of this background has to be explored more thoroughly in the future. 

In the close vicinity of the spike, one notices that in the upward direction 
the red field distribution has a long tail. For the green colour, one observes a 
flip in the orientation at a distance A from the center of about 3 lattice units. 
Furthermore, the blue fields flip sign when moving right. Thus, a possible 
instanton wave function is strongly distorted for A > 3a, if it is realized in the 
regular gaugeQ 

It has to be tested whether an anti-instanton realized in the singular gauge 
provides a better description. Consequently, in fig. |9| I show the fields of a broad 
artificial instanton with p = 15 a in singular gauge, subject to the same F.f. The 
field lines are approximately parallel straight lines which flip orientation at a 
distance of A = 5a. This comes as an interplay of the singular gauge and the 
cut-oflQ. This flip is also noticable in the previous figure for some colours and 
some directions. A scan through all configurations and directions reveals that 
this phenomenon is frequent. It is one of the indications that we are dealing 
with singular gauges and less likely with pointlike instantons. On the other 
hand, the long tail of the field with red colour observed in the previous figure, 
is also not reproduced by this artificial anti-instanton. 

In fig. H I display field lines. Starting from a local maximum of S(x), electric 
field lines for all colours have been drawn, where the field strengths off the lattice 
sites have been determined by interpolation. When rotating this picture with 
the help of the VRML-browser, on the one hand, one clearly checks for the 
orthogonality of the field lines close to the point of intersection, but, on the 
other hand, one realizes that after 4 cones (two lattice units) the field lines 

7 Significant fluctuations at a distance of less than 3 lattice units are not expected to occur, since the 
cut-off suppresses Fourier strongly amplitudes with k 2 > 1/a 2 . 
8 For a stronger cut-off, the field lines become more parallel. 
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depart from the original direction. Again, this means that the wave function of 
the object no longer agrees with that of an instanton. A plot for the magnetic 
field lines shows exactly the same trends, consistent with selfduality. 

Proceeding to larger values of the cut-off, i.e. A 2 = 4 a 2 , I show, in fig. [10], the 
electric fields of red colour in a plane through the peak. Now, deviations from the 
instanton configuration (parallel field lines) occur at a distance A of 10 a, and 
at a distance A = 15 a one can observe a few other, less pronounced maxima. 
This characteristic distance has to be interpreted in connection with the cut-off 
in momentum space. If we take the condition kA = n for the occurence of a 
significant deviation and require X 2 k 2 < 0.5, we obtain A > 8.9 a, in agreement 
with the visual impression. In other words, the wave function is disturbed as 
soon as allowed by the cut-off. 

Next, in figs. |TT] to ^ the gauge field distributions in the same plane as 
previously are shown, again for A 2 = 4 a 2 . Now, the peaks in the action and in 
P a (x) are clearly correlated in space to a multi-dimensional zero in the gauge 
fields, i.e. a zero for all directions and for all colours. This zero has evolved in 
quality when coming from smaller values of the cut-off, where it is obscured by 
the singular nature of the unsmoothed gauge fields. 

We note that the zero is coupled to regions of large gauge fields. However, 
there are several other excitations of similar magnitude visible. Furthermore, 
the environment of the instanton does not seem to be the only region where rapid 
variations of the gauge fields occur. Among others, one notices the occurence 
of several vortex-like objects, some of those within a distance A < 15 a. The 
correlation of the action density with these objects remains to be investigated. 

Whereas the spikes are phenomena with a well defined density and, apart 
from the radial dependence, with clear properties of instantons, the maxima 
of the "background" are more difficult to interpret. From fig. |] we conclude 
that they are much more numerous than the spikes, and that they do not obey 
scaling of the radius, by the same arguments as given for the spikes. To isolate 
possible candidates of instantons with larger radii will require different filtering 
techniques as have been used up to now. 



5 Cooling vs. Fourier Filtering 

In this section, I study the connection between the peaks emerging after Fourier 
filtering, and peaks emerging after cooling. For both smoothing techniques, I 
start from the same configuration in the Landau gauge. There are two pos- 
sibilities for a comparison. Either, one observes the gauge fields after cooling 
around those lattice points which are selected in the search for maxima in F.f., 
or one searches for maxima in the cooled configuration independently and plots 
the environment around those positions. 

Before discussing the outcome in detail, the cooling technique has to be 
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specified. There one replaces the link variables iteratively by a weighted sum of 
the link variables U^(x) and the staples Tl^x), 

%{x) = V((l - c)U,(x) + |e„(x)), (27) 

where V stands for projection onto SU(2)-matrices. The links are updated 
following an arbitrary path through the lattice. Based on comparison with 
renormalization group mapping, a value of c ~ 0.5 is recommended Q. 

It is obvious that this method, after sufficiently many iterations, will drive 
the lattice configuration to a superposition of local action maxima^] which are de- 
termined by the local boundary conditions (generated by concurrent maxima), 
thereby minimizing the global action. These maxima may be instantons. The 
density of instantons and anti-instantons, generated in this way, is a strongly 
decreasing function of the number of cooling steps [Q, [7] , whereas the difference 
between the number of instantons and anti-instantons depends not so strongly 
on cooling. This allows to determine the phenomenologically interesting quan- 
tity < Q 2 > within reasonable errors. 

To start a comparison, I refer to the global plot of the pseudoscalar density 
P a (x) of the configuration used in figs. [I] and || on the large lattice at (3 — 2.85. 
There, the positions of the highest peaks, obtained after F.f., had been selected, 
and the cut-off was A 2 = 0.3 a 2 and A 2 = 1.0 a 2 resp. In fig. [14| the same config- 



uration is shown, after 30 cooling sweeps with c = 1 (see eq. \27\). The cooling 
started from Landau gaugeQ, and the maxima were selected after cooling. The 
scale of the cones has been increased by a factor 2.5 as compared to fig. |2|. 

Obviously, the overall characters of the plots, obtained by F.f. and by cool- 
ing, are quite different. In the latter case, the peaks are much broader and 
have smaller amplitudes at their centers. The alignment along the diagonal in 
colour space is much more pronounced, and there is only one prominent peak 
as compared to the six objects visible in fig. (remember the different scale in 
the two figures!). 

From these plots, a direct comparison between peaks, eventually located at 
the same positions, is not easy. In order to compare corresponding peaks in 
the two cases, it is necessary to specify the positions of the peaks obtained by 
F.f. and to plot P a (x), after cooling, at these locations. The outcome of this 
procedure is the following: One again finds maxima around the positions of the 
peaks, apart from a shift of the location of the maxima by a few lattice units. 
The spatial extension of the peaks is much larger, their amplitudes are reduced 
by a factor rs 100, and, most importantly, the sign of P a (x) has changed. The 
origin of this behaviour is clearly demonstrated in fig. [15], where the gauge fields 



3 The term maxima refers to local spatial maxima, and cooling brings the configuration to global action 



minima 
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This is irrelevant for the discussion of the variable P a (x) 1 but not for the comparison of the gauge 
fields. The gauge independence of the action and of the pseudoscalar density will be exploited in the next 
section. 
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with red colour are displayed^ around the position of the peaks found in F.f. 
Obviously, by cooling the configuration evolves towards a gauge singularity, 
where the gauge fields are maximal in absolute size and vary most rapidly 
(changing sign) at the position of the peaks. 

It is instructive to monitor the influence of cooling and F.f. by an animated 
presentationP]. After each cooling step and after each increase of A 2 by 5\ 2 = 
0.05 a 2 , the corresponding gauge fields and P a (x) around a selected peak are 
written to a file. The VRML language allows to convert the cooling/filtering 
steps into a time evolution. Depending on the number of sites in the plot, 
the change of the gauge fields in time is reasonably smooth, due to the linear 
interpolation offered by VRML. One observes that the gauge fields far away (by 
3 lattice units or more) from the site of the peak decrease, in both schemes, 
more or less in the same way. However, the gauge fields which are located 
closely to the site remain in size essentially unchanged under cooling, with a 
few link variables rotating through colour space. This leads to the singularity 
in the gauge fields. As an effect of the minimization, the action and P a (x) 
are quite small. This comes about via the cancellation of the linear and the 
quadratic terms in the field strengths^]. In F.f., on the other hand, the decrease 
in amplitude is uniform, and the gauge fields develop their smooth zero at the 
point of the spike. Now, the action and P a (x) are not small, since the typical 
cancellation mentioned before is not operative. 

It has already been mentioned that far away from the spike, i.e. for A > 3 a, 
the gauge fields are similar in both schemes (in detail, the fields are somewhat 
more rough after cooling than in F.f., and their absolute magnitude depends on 
the amount of cooling and, in the case of F.f., on A 2 ). Thus, one can say that the 
two schemes lead to a different interpolation between regions of blocked gauge 
fields. In F.f., this connection is a smooth one with a zero at the peak position, 
leading to an approximate instanton (or anti-instanton) wave function with 
small radius p in the regular gauge. After cooling, the connection is achieved by 
a singularity, leading to the wave function of a broad anti-instanton (instanton) 
in the singular gauge. Probably, it is very difficult to decide which scenario is 
the "correct" one. In principle, the behaviour of the gauge fields at A « la 
and at A > p should provide an unambiguous tool for distinguishing between 
the two cases. However, in the first region, lattice artifacts prevail, and in the 
second region the wave functions are distorted heavily after a couple of lattice 
units, as can be seen in fig. 

Finally I note that the outcome of cooling is influenced by high momentum 
amplitudes. If one applies, before cooling, a F.f. with a small value of the cut- 

n It has to be emphazised that the colour in this representation indicates the position in the fourth 
dimension. 

12 Unfortunately, it seems to be impossible to demonstrate this in a publication. 

13 This can be easily demonstrated by retaining only the linear term in the fields, in which case one 
observes the same spikes as in F.f. 
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off, A 2 = 0.2 a 2 , the configuration evolves to a smooth field distribution when 
being cooled. 

In summary, both cooling and F.f., when applied to a configuration in Lan- 
dau gauge, lead to a set of points (with a density of 1.4 points per fm A ) where 
the gauge fields in the environment are connected by a rapid flip of signs for 
all amplitudes. These points signal the impossibility to find a gauge which is 
smooth, after the removal of perturbative fluctuations. When one does the lat- 
ter with the help of cooling, the connection between regions of different sign is 
achieved by a singularity in the gauge fields, and this allows for a smooth and 
minimal distribution of the action. After F.f., the gauge fields are interpolated 
smoothly (details depend on the cut-off), and this leads to a spike in the action. 
The latter fact is due to the inadequacy of linear modifications at the singular 
point, such that the subtle cancellation does not work. Anyhow, for both meth- 
ods, visualization is a good tool to study the evolution of variables as a function 
of the number of cooling steps or of the amount of Fourier filtering. 

6 Gauging Ambiguity 

Gauging is an essential tool in the present investigation, and a natural reser- 
vation is based on the corresponding ambiguities. The latter are given by the 
type of gauge (Landau gauge vs. Coulomb, axial or temporal gauge) and by the 
Gribov ambiguities Jl| within a particular gauging scheme. 

6.1 Gribov Ambiguities 

In order to study the influence of the Gribov ambiguities, five configurations at 
/3 — 2.7 have been 12 times subject to the following procedure: Before trans- 
forming to the Landau gauge, a random gauge transformation has been applied 
to the configuration. In about 20% of these 12 trials, the ensuing Landau gauge 
fixing algorithm ended up in one specific gauge. For each gauge, local maxima 
of the action density have been determined in the same way as described in 
section £|, and the positions of the 15 strongest maxima have been compared. A 
shift of 2 lattice units in every direction was allowed for considering the positions 
as identical. For all 5 configurations, there is at least one spike which occurs at 
the same position in all of the 12 gauges. For about 20% of the maxima, one 
gauge missed one of the spikes, and for 2 misses, there is a similar probability. 

It has to be concluded, that the positions of the spikes are not completely 
unambiguous properties of the gauge orbit of a given lattice configuration. But 
they are also not fully a gauge artefact either. Apparently, the configuration 
has, at certain positions, a distortion which acts as an obstruction against the 
interpolation by a smooth gauge. The gauging algorithm either may find a 
gauge field realization which leads to strong gradients at these positions, or a 
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smooth one, shifting the problem to another region of space-time. 

It is interesting to note that in the Coulomb gauge the existence of gauge 
singularities is tantamount to the occurrence of Gribov ambiguities fl8f . Since 
we can choose the coupling constant go small enough such that no spikes are 
present in our configurations, it will be possible to test the connection between 
the phenomena numerically, albeit in Landau gauge. 

6.2 Spike Positions in F.f. and in Cooling 

The gauge independence of the action and of the pseudoscalar density under 
cooling offers another possibility to study the gauge dependence of the F.f. 
method. For several configurations, I have applied 15 cooling steps with c = 
0.5, and afterwards the pseudoscalar density P a (x) has been plotted on planes 
running through spikes found under F.f. These plots can be compared directly 
with those obtained by F.f., e.g. with a cut-off A 2 = 1.0 a 2 . In almost all cases it 
turns out that at the position of the spike there exists a maximum of | P a (x) | in 
the cooled configuration. Furthermore, a strong correlation exists also between 
most of the secondary peaks in F.f. and in cooling. The spatial extension of 
these peaks as well as the orientation of P a (x) in colour space differs in the two 
cases. This is not surprising. After all, there is a strong dependence of these 
quantities both on the number of cooling steps and on the cut-off used in F.f., 
and different quantities are being optimized in the course of the distortion of 
the lattice configuration. 

One can conclude that the gauge dependent peaks obtained by F.f. show up 
at essentially the same positions as those emerging from the gauge independent 
cooling method. This makes it quite unlikely that in F.f. we are dealing with 
gauge artefacts. However, the physical properties of the peaks are necessarily 
different in the two methods. 



7 Comparison with QED 

How can the vacuum of SU(2) lattice gauge theory be distinguished from the 
free field case, i.e. from non-compact QED? The latter will simply be defined 
by a superposition of plane waves on the lattice. After creating such configura- 
tions, they will be subject to the same Fourier filtering as the SU(2)-vacuum, 
without encountering the problems of linearization, of course. For QED, no 
gauge singularities are expected and none are found. Thus, both vacua differ 
in a fundamental aspect. However, if a lattice plane is selected which does not 
pass through a spike, it is impossible to differentiate between both theories by 
inspecting plots of the action or of the gauge fields. An example is given by 
a comparison of fig. [16| with fig. [17], for a cut-off A 2 = 2 a 2 and a lattice size of 



48 3 x 64. For SU(2), (3= 2.85. The author has not been able to pin down a 
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significant difference in the size of the gauge field clusters, which survive the 
smoothing procedure. 

In order to make the comparison more quantitative, the correlation of gauge 
fields around the various maxima is studied. For a fixed colour a, the gauge 
fields around many (numerically determined) maxima are averaged within a 
distance of 2 lattice units. This average will be denoted by A (x). The same 
average is taken around a center shifted in the direction of the local field by a 
distance R. The scalar product of the two averages is taken, and an effective 
longitudinal mass mi(R) is determined by 

mi (R) = -\nJ2Al(x)Al(x + e A R)/j:(A;( x ))^ ( 28 ) 

where e& is a unit vector into the direction of A a ^(x). Because we move around 
a maximum of the square of the gauge fields, the slope at small values of R is 
small. In a similar way, transverse correlations are defined for the gauge fields. 
In fig. the transverse and longitudinal masses both for QED and for SU(2) 
are compared. The lattice size is 32 3 x 64, and for SU(2) I take (3= 2.7. In 
agreement with expectation, the masses in SU(2) are significantly smaller than 
in the case of QED. The reason for this phenomenon lies in the behaviour of 
the gluon propagator in Landau gauge PJJ, E2| which increases more strongly 
towards small momenta than the free propagator Go(q 2 ) = 1/q 2 . Transformed 
to x-space, this implies a softer decrease of the correlation in SU(2) as compared 
to QED. Although the differences in the correlation are significant, they are too 
small as to be visible in a single configuration. 

Not so clear is the behaviour of the correlations at large R. Especially in 
the transverse case, the correlation in QED seems to develop a zero around 
R = 7a, which is, of course, cut-off dependent. In SU(2), no zero is found. An 
explanation of this phenomenon is outside the scope of this paper. 

Based on the picture of stochastically oriented domains of gauge fields [ITU 



20| , I had expected to obtain a rapidly vanishing correlation in the non-abelian 
case, i.e. an increase of the effective masses as functions of R. This has not yet 
been observed. Presumably, the distances R used here are too small as to see 
this effect. 



8 Conclusions 

There is no doubt that computer (PC-) technique has matured sufficiently such 
that an instructive visualization of large sublattices in SU(2) lattice gauge theory 
is feasible. Still, the huge number of variables imposes some restrictions in the 
presentation of lattice variables, but these are not so strong as to render the 
information obtained by visualization meaningless. This is demonstrated by the 
observations described in the previous sections, which were new at least to the 
author. In recapitulation, I found 
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1. the existence -in Landau gauge- of objects with a rapid variation of gauge 
fields which persist through Fourier filtering up to large values of the cut- 
off. Under this technique, they show up as narrow peaks in the (gauge non- 
invariant) action and pseudoscalar density. The density of these objects is 
1.4/m~ 4 , it scales according to the string tension, and it agrees with the 
density of instantons found by a special cooling method 0. 

2. The properties of these objects, e.g. the orientation of electric and magnetic 
field lines within the peaks, both in x-space and in colour space, are those 
of instantons. 

3. The peaks are well separated from a background which does not resemble 
the characteristics of instantons in any obvious way, at least according to 
first impression. More work is necessary here. 

4. The strong concentration of field strenghts in x-space is compatible either 
with the existence of very narrow (p w la) instantons in regular gauge or 
with the existence of broad instantons (p > 10 a) in singular gauge. The 
latter case is favoured by the fall-off of the gauge fields as a function of the 
distance from the center. 

5. The quality of the wave function of the instantons, as measured by the 
field strengths, is strongly corrupted for distances from the center, which 
are smaller than the radius < p >, quoted in the literature. 

6. The different response of the gauge field configuration and of other variables 
either to Fourier filtering or to cooling can be demonstrated impressively 
by animation. 

7. The positions of the spikes are partly gauge dependent, i.e. some differ 
for different Gribov copies. However, all spikes are associated with action 
maxima which appear after cooling, some of which have quite small values 
of the action close to the center. 

8. The similarity of QED and SU(2) in their respective gauge field distri- 
butions, apart from the existence of the gauge singularities, can be made 
visible. 

It is to be expected that with more experience additional insight into other 
structures of the vacuum will be obtained. 
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Figure Captions 

The figures (mostly coloured) are in the format .png, except the two black-white 
ones, which are in the .ps format. 



Figure 1: Large lattice (48 3 X 64, (5 = 2.85): Pseudoscalar density maxima, 
projected into 3 dimensions. Colour indicates position of maximum in fourth 
dimension, moving from red (z = 1) to green (z = 48). The cut-off is A 2 = 0.3a 2 . 
The direction of cones indicates the direction in colour space. 



Figure 2: Large lattice: Pseudoscalar density maxima, projected into 3 di- 
mensions. Colour indicates position of maximum in fourth dimension, moving 
from red (z = 1) to green (z = 48). The cut-off is A 2 = 1.0a 2 . The direction of 
cones indicates the direction in colour space. 



Figure 3: Electric field lines passing 
colours are those of SU(2). The scale 
distance between 2 consecutive cones is 
1.0. 



through a maximum of action. The 
is such that close to the center the 
half a lattice unit. Large lattice, A 2 = 



Figure 4: Histograms for summed Wilson action E(A 2 ,i?) (see eq. ^) for f3 
= 2.85, size = 48 3 x 64 (full line), and (3 = 2.70, size = 32 3 x 64 (broken line). 
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Figure 5: Small lattice (32 3 x 64, (3 = 2.70): Pseudoscalar density maxima, 
projected into 3 dimensions. Colour indicates position of maximum in fourth 
dimension. The cut-off is A 2 = 1.0 a 2 . 



Figure 6: Large lattice: Pseudoscalar density P a (x) in a lattice plane running 
through a spike. The cut-off is A 2 = 0.15a 2 . The direction of the cones is given 
by the orientation in colour space. 



Figure 7: Large lattice: Gauge fields of red colour in a lattice plane running 
through the spike shown in the previous figure. The cut-off is A 2 = 0.15 a 2 . 
The direction of the cones is given by 3 space-time components of the gauge 
fields. 



Figure 8: Large lattice: Electric fields around a spike. The cut-off is A 2 = 
0.5 a 2 . The colours are those of SU(2). 



Figure 9: Electric fields around an artificial instanton with radius p = 15 a 
in singular gauge. The cut-off is A 2 = 0.5 a 2 . The colours are those of SU(2). 



Figure 10: Large lattice: Electric field strengths of red colour in a lattice 
plane running through the spike shown in fig. [(]. The cut-off is A 2 = 4.0a 2 . The 
direction of the cones is given by the spatial orientation of the fields. 



Figure 11: Large lattice: Gauge fields of red colour in a lattice plane running 
through the spike shown in fig. |6|. The cut-off is A 2 = 4.0a 2 . The direction 
of the cones is given by three space-time directions of the gauge fields. The 
magnitude of the pseudoscalar density is indicated by golden circles. 
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Figure 12: Large lattice: Gauge fields of green colour in a lattice plane running 
through the spike shown in fig. |5|. The cut-off is A 2 = 4.0a 2 . The direction 
of the cones is given by three space-time directions of the gauge fields. The 
magnitude of the pseudoscalar density is indicated by golden circles. 



Figure 13: Large lattice: Gauge fields of blue colour in a lattice plane running 
through the spike shown in fig. ||]. The cut-off is A 2 = 4.0a 2 . The direction 
of the cones is given by three space-time directions of the gauge fields. The 
magnitude of the pseudoscalar density is indicated by golden circles. 



Figure 14: Large lattice: Pseudoscalar density, projected into 3 dimensions. 
Colour indicates position in 4th dimension. 30 cooling steps have been ex- 
ecuted with c = 0. The positions of the peaks are found numerically after 
cooling. 



Figure 15: Large lattice: Gauge fields of one SU(2) colour, at positions of 
spikes which are determined by Fourier filtering. Positions are projected into 
3 dimensions, and colour indicates 4th dimension. The configuration has been 
cooled by 30 steps with c = 0., after transforming it to Landau gauge. 



Figure 16: Large lattice: Gauge fields (3 components only) in a plane for 
SU(2), avoiding spikes. The cut-off is A 2 = 2.0a 2 . 



Figure 17: Gauge fields in a plane for non-compact QED. The cut-off is 
A 2 = 2.0a 2 , and the lattice size is 48 3 x 64. 
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Figure 18: Effective masses for longitudimal and transverse correlations of gauge 
fields. From top to bottom: Tranverse QED (open circles), transverse SU(2) (open 
squares), longitudinal QED (open diamonds), longitudinal SU(2) (full circles). The 
lattice size is 32 3 x 64, the Fourier filtering parameter A 2 = 1.0 a 2 . For SU(2), (3 = 
2.7. 
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This figure "il.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/00 1 1 046v 1 



This figure "i2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/00 1 1 046v 1 



This figure "i3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/00 1 1 046v 1 
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Figure 4: Histograms for summed Wilson action E(A 2 , R) (see eq. 23) for (3 = 
2.85, size = 48 3 x 64 (full line), and (3 = 2.70, size = 32 3 x 64 (broken line). 
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Figure 18: Effective masses for longitudimal and transverse correlations of gauge 
fields. From top to bottom: Tranverse QED (open circles), transverse SU(2) (open 
squares), longitudinal QED (open diamonds), longitudinal SU(2) (full circles). The 
lattice size is 32 3 X 64, the Fourier filtering parameter A 2 = 1.0 a 2 . For SU(2), (3 = 2.7. 



